State-dependent lattices for quantum computing with alkaline-earth-metal atoms 
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Recent experimental progress with Alkaline-Earth atoms has opened the door to quantum com- 
puting schemes in which qubits are encoded in long-lived nuclear spin states, and the metastable 
electronic states of these species are used for manipulation and readout of the qubits. Here we dis- 
cuss a variant of these schemes, in which gate operations are performed in nuclear- spin- dependent 
optical lattices, formed by near-resonant coupling to the metastable excited state. This provides 
an alternative to a previous scheme [A. J. Daley, M. M. Boyd, J. Ye, and P. Zoller, Phys. Rev. 
Lett 101, 170504 (2008)], which involved independent lattices for different electronic states. As in 
the previous case, we show how existing ideas for quantum computing with Alkali atoms such as 
entanglement via controlled collisions can be freed from important technical restrictions. We also 
provide additional details on the use of collisional losses from metastable states to perform gate 
operations via a lossy blockade mechanism. 



I. INTRODUCTION 

There has been a lot of recent experimental progress 
in cooling and manipulating alkaline-earth and alkaline- 
earth-like atoms in the laboratory, especially in the con- 
text of optical clocks with Strontium Atoms [IH5], and 
the production of Bose-Einstein condensates and degen- 
erate Fermi gases of Ytterbium [SHE]) Calcium [3] and 
Strontium [TTJMl3| . The control that has been developed 
over these atoms makes them an extremely interesting 
candidate for the implementation of quantum informa- 
tion processing [TH - fTT] . This is especially true in light of 
the laser stability achieved in optical clock experiments 
[S] , which is reminiscent of the development path to- 
wards quantum computing taken in the case of trapped 
ions QUO!]. 

The key new feature of alkaline earth atoms in com- 
parison with alkali atoms is the singlet-triplet metastable 
transition, with the 1 So - 3 Po transition being used as 
the clock transition (see Fig. [l^i). In particular, for 87 Sr, 
the 3 Po manifold has a measured lifetime of r ~ 30s, 
and the 3 P2 levels have even longer predicted lifetimes. 
In addition, for species with non-zero nuclear spin, this 
spin can be decoupled from the electronic state on the 
clock transition [TAHITI 12U] , especially in the presence of 
a large magnetic field. The use of this nuclear spin for 
storage of quantum information would then be ideal, as 
the nuclear spin is much less sensitive to magnetic fields 
than electron spins, and thus much less susceptible to 
decoherence from magnetic field fluctuations than qubits 
stored on electronic states. This has lead to a series of 
proposals [HI Q~51 HI] in which the electronic state is used 
for access to and manipulation of the qubit |22j , and the 
nuclear spin state is used for qubit storage. 

In previous work [14] we developed a scheme for quan- 
tum computing with alkaline earth atoms that was based 
on electronic- state- dependent lattices in which indepen- 




FIG. 1: Level structure for Alkaline-earth-like atoms, (a) 
These atoms possess a singlet-triplet transition with long- 
lived metastable 3 Po and 3 P2 levels, (b) Adiabatic dressed 
potentials can be created by a resonant coupling on the clock 
transition with a sinusoidally varying Rabi frequency Qi(x), 
producing dressed states of the 1 So and 3 Po levels, (c) The 
differential Zeeman shift for different nuclear spin levels can 
be used to produce nuclear spin-dependent lattices, by driv- 
ing the transition between the 1 So and 3 Po levels resonantly 
at the different resonant frequencies for different nuclear spin 
states in a magnetic field. Here we show some of the mi states 
for an atom with nuclear spin I = 9/2, such as 87 Sr. 



dent control over lattices for the metastable excited 3 Po 
and ground 1 So levels is obtained by using light of differ- 
ent wavelengths. This is made possible by the fact that 
these levels are optically separated, providing very dif- 
ferent AC polarisabilities for the states as a function of 
the wavelength. We showed how these two independent 
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lattices could be used as a storage lattice for qubits en- 
coded on the nuclear spin state, and a transport lattice 
to manipulate the qubits and perform gate operations 
[13] . A key theme in this context is that many schemes 
and concepts developed for alkali atoms, including cer- 
tain techniques that have already been demonstrated in 
proof-of-principle alkali experiments simply work quan- 
titatively better for alkaline earth atoms, where they are 
freed from important technical restrictions. In this sense, 
alkaline earth atoms represent an extremely important 
technological advance in various settings. 

Here we present an alternative scheme to this previ- 
ous proposal, in which we show that near-resonant cou- 
pling on the clock transition can produce frequency selec- 
tive nuclear-spin-dependent lattices (see Fig.[l]D,c). As in 
the case of electronic-state-dependent lattices, this makes 
it possible to produce state-dependent lattices without 
the restriction of having to tune couplings between fine- 
structure states [53J |23J, which can lead to large heat- 
ing and decoherence from spontaneous emissions in the 
case of alkali atoms. In addition, spin-dependent lat- 
tices made in this way can be easily generated so that 
motion of the two potentials is independent in 2D. Be- 
low we discuss this implementation in detail, develop- 
ing a complete proposal for quantum computing with 
alkaline-earth(-like) atoms, including methods for pro- 
duction of a quantum register and for performing single- 
qubit operations. Qubit readout with individual address- 
ing can be performed in a similar manner to the case of 
electronic-state-dependent lattices [T3J, using magnetic 
gradient fields to shift the energy of states in the 3 P2 
manifold. In nuclear-spin-dependent lattices, the large 
two-body loss rates from metastable 3 P2 levels can also 
be used to perform two-qubit gate via a lossy blockade 
mechanism as an alternative to the implementation of 
gates via controlled collisions [23]. This was originally 
discussed for electron-state-dependent lattices [14] . and 
we provide further details of this mechanism below. 

This alternative scheme with nuclear-spin-dependent 
lattices has the advantage of not requiring additional 
lasers to trap the 1 S and 3 Po lattices independently, 
and can be performed with a strong laser on the clock 
transition. At the same time, this method is somewhat 
sensitive to collisional losses when two atoms in the 3 Po 
manifold collide, as the near-resonant lattices will always 
produce admixtures of this state. However, this is only 
an issue during the short times involved in gate oper- 
ations, and is strongly suppressed in a realistic setup 
where atoms are also AC-Stark shifted, as discussed be- 
low. Nuclear-spin-dependent lattices would also have im- 
mediate applications in quantum simulation [25 -29 . In 
particular, the degeneracy in models with SU(N) symme- 
try [2"5Tf2"T] (which can be studied using alkaline-earth-like 
atoms by making use of the symmetry for interactions of 
atoms in different nuclear spin levels) could be deliber- 
ately broken and restored by applying these nuclear-spin- 
dependent potentials. 

The rest of this article is organised as follows: We first 



discuss the formation of near-resonant spin-dependent 
optical lattices in more detail in Sec II, together with 
preparation of atomic registers in these lattices. In 
Sec. Ill we then discuss means for readout of individual 
qubits, and in Sec IV we treat ideas for gate schemes to 
entangle two qubits, including making use of lossy block- 
ade mechanisms. In Sec. V we present a summary and 
outlook. 



II. SPIN-DEPENDENT ADIABATIC 
POTENTIALS 



In the following, we discuss a quantum register formed 
by one atom trapped every site of a deep optical lattice, 
where tunnelling of atoms between sites can be neglected 
on the timescale of the experiment. As discussed above, 
we identify hyperfine states with two chosen nuclear spin 
states, and we would like to create spin-dependent poten- 
tials in order to move these qubit states independently. 
We will make use of these in the two-qubit gate opera- 
tions that we discuss in Sec. IV. 

A novel method of forming optical lattices for alkaline 
earth atoms is to make use of a near-resonant optical 
coupling directly on the clock transition, which will pro- 
duce adiabatic dressed potentials [S3]. In the case that 
the coupling field is a standing wave, the Rabi frequency, 
and thus the final dressed potential, will be sinusoidally 
varying, providing an optical lattice for dressed states 
that are superposition of states in the 1 So and 3 Po levels 
(as shown in Fig. [ljs). In a large magnetic field, there 
is a differential Zeeman shift AEz between the 1 So and 
3 Po states (109 Hz/G for 87 Sr [5]), meaning that a direct 
coupling preserving the nuclear spin (with 7r-polarised 
light) will be resonant at substantially different frequen- 
cies for different nuclear spin states (see Fig. [TJ:). We 
can then drive each transition independently with Rabi 
frequencies Oq = as shown in Fig. [T]d. Provided 
that the shift AEz 3> £l±, we will then obtain indepen- 
dent two-level systems for each mj state for which we 
apply the appropriate coupling frequency. For example, 
if we choose Slo ~ 100 kHz, then for 87 Sr, we would like 
to apply a field > 1000 G in order to obtain shifts be- 
tween neighbouring states > 100 kHz (however, states 
separated further in mj could also be used to reduce the 
required field - see below). In this way, we can choose, 
e.g., two m; states as our two qubit states, |0) and |1), 
and create independent potentials for these two states. 
At the same time, because the frequency differences be- 
tween lattices for different mi states will of the order of 
1 MHz, the lattice laser wavenumber ki is approximately 
the same for the two species - in fact the resulting lattice 
potentials will overlap for the order of millions of periods. 
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A. Dressed potentials for a two-level system 

We will now discuss the form of the dressed potentials 
for a single nuclear spin state, identified with qubit state 
i = or i = 1, and discuss the case where we have 
multiple nuclear spin states below. We can first write 
the Hamiltonian for a two-state atom, with states \g, i) = 
pSo, mi = i) and |e, i) = | 3 Po, i) as (h = 1) 

H = H M + H Q , 

where Hm = p 2 /2m is the kinetic energy, and 

H = -Si\e,i){e,i\ + (fii(x)/2)|e,t)<$,t| + h.c. 

describes the near-resonant coupling field with £li (x) and 
Si the Rabi frequency and detuning respectively. 

Generation of adiabatic dressed potentials is then 
based on the validity of a Born-Oppenheimer-type as- 
sumption, where we assume that the kinetic energy of 
the atoms is small on a scale given by the separation 
of the resulting adiabatic potentials. The wavefunction 
|$(£)) of a single atom satisfies the Schrodingcr equation, 

If we omit the kinetic energy term from the Hamiltonian, 
we obtain an equation for adiabatic eigenstates, \^±(t)), 

H (x)\* ± )=V ± (x)\* ± ), (2) 

Note here that as i?o( x ) is time-independent, there are 
only two such eigenstates \^f±). If we consider the ID 
case, and set ili(x) — ttisin(kix + </>), representing the 
field of a standing wave (with ki the laser wavenum- 
ber and <j) a phase), we find the adiabatic potentials 
V ± {x) = (Si ± y/S^ + Qi(x) 2 )/2. These are shown 
schematically in Fig. [Tb). The complete wavefunction 
can then be expanded in a basis of these adiabatic eigen- 
states, which play the role of Born-Oppcnheimer channel 
functions, 

m)) = c + (x,t)\* + )+c_(x,t)\*-), (3) 
resulting in the equation 

ih~c±(x,t) = [H M + V ± (x)]c ± (x,t) + H±c T (x,t), (4) 

where = fif±(x,t)\HM\&^(x,t)) gives the non- 

adiabatic couplings between the dressed states due to 
the motion of the atom. Provided these latter terms are 
small, the two equations decouple and the atoms remain 
in a single dressed state. In our case, there will be no 
non-adiabatic loss of atoms in this sense, provided that 
they are loaded into low energy states of the lower of 
the two adiabatic potentials (V~). If atoms are loaded 
into the higher energy dressed potential, loss of atoms 
into the continuum states of the V~ potential can occur, 
however this is exponentially suppressed as the separa- 
tion between adiabatic potentials is increased, with the 
loss rate T[ ~ T exp(— <5j/w) [20], where T is a prefactor 
that we do not compute in detail here, and u> is the trap 
frequency in an individual site of the dressed potential. 



B. Dressed potentials for independent states 

When these dressed potentials are created indepen- 
dently for different mi levels, the result is that we can 
form two independent but almost identical potentials. 
These can, e.g., be shifted with respect to each other 
in a 2D plane using interferometrically stable methods, 
e.g., by adding path length to an interferometer arm in 
which the light is frequency-shifted in order to produce 
one of the trapping frequencies. As discussed above, this 
means of creating spin-dependent lattices has substantial 
advantages over spin-dependent lattices for Alkali atoms, 
where the lifetime is limited by the need to tune the lat- 
tice beams to a frequency in the middle of the fine struc- 
ture splitting. Here, the lifetime will be controlled by the 
lifetime of the 3 Po level (which is many seconds), or by 
off-resonant couplings to shorter lived states (but these 
will typically be many tens of nanometers detuned). In 
the presence of a second frequency (e.g., due to the laser 
creating the lattice for the second internal state), atoms 
can also be lost from the lower adiabatic potential, essen- 
tially being coupled out of the lattice into the continuum. 
However, due to the large momenta in the resulting state, 
this rate is suppressed exponentially in the ratio of the 
separation between manifolds in a Floquet basis and the 
trapping frequency in the lattice, T; ~ To exp(— Wdifi/w) 
[5U] , where r is a prefactor [2U] Wdig is the frequency dif- 
ference of the lattices for the two qubit states. If we oper- 
ate in a field ~ 5000 G, then Wdifi ~ 27r x 550kHz, and if 
we choose the Rabi Frequency 51 ~ 120 kHz, then u ~ 15 
kHz. We can also reduce the required field strength by 
choosing mi levels that are further separated (in 87 Sr we 
can reduce the required field strength by a factor of 9 by 
choosing mi = —9/2 and m; = +9/2 as the two trapped 
states) . 

C. Combining resonant and off-resonant potentials 

In practice, strong coupling at intensity / on the clock 
transition will also give rise to off-resonant AC-Stark 
shifts AE AC and &E 9 AC of the states \e,i) and \g, i) 
from coupling to other manifolds (e.g., 1 Pi and 3 Si) 
in addition to the resonant couplings between the two 
levels. These must be added to the Hamiltonian, as 
H = H M + H + AE AC \e,i){e,i\ + AE AC \g,i){g,i\. As 
AE A q oc I and fl cx y/l, the off-resonant contributions 
will become more important as the intensity of the ap- 
plied field becomes larger. For 87 Sr, the shifts from the 
AC-Stark shift become of the same order as the AC-Stark 
splitting due to resonant coupling at relatively high fields, 
with / - 50kW/cm 2 [20]). At higher fields, the po- 
tentials will be modified by these shifts, but can still 
be made spin-dependent if the detunings and Rabi fre- 
quencies of the lattice beams are chosen carefully. This 
is illustrated in Fig. [2j where we show the lower adia- 
batic potential for each of the two nuclear spin states for 
a selection of different phases <f> between the potentials. 
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FIG. 2: a) Lower dressed potentials for the |0) qubit state 
(solid lines) and the |1) qubit state (dashed lines) formed by 
combining resonant and off-resonant contributions. These are 
plotted as a function of position for varying phase offsets c/> 
between the coupling fields for the |0) and |1) states. For 87 Sr 
at the wavelength of the clock transition, we obtain AE% C w 
3 AE 3 AC , and we choose the intensity so that Q — 4AE AC (ca. 
10 kW/cm 2 ). Here, 5 — — 3J7/4. b) Projection on the excited 
state of the dressed state corresponding to the lower dressed 
potential, plotted for the same values of Q, AE AC , and AE AC 
as in (a) as a function of c/>. The different lines correspond to 
varying S, from top to bottom, S — — $7/2, — 3f7 /4, — Q, and 
-5J7/4. 



We see that at relative phase = the potentials for 
different nuclear spin states are identical, and are given 
by a combination of the resonant and off-resonant con- 
tributions. At phase <j) = it/ 2, however, the off-resonant 
contributions from the two coupling frequencies, which 
are independent of the nuclear spin state, become spa- 
tially homogeneous due to the addition of the two spa- 
tially shifted contributions. At this point the sinusoidal 
form of the lattice potentials is due solely to the resonant 
contribution. It can be seen that the lattice, also in be- 
tween, will be modified in such a way that the atoms will 
be transported through the lattice spin-dependently. 



We also note that two important characteristics re- 
lating to the shape of the lattice and the form of the 
dressed states changes as a function of 0. Firstly, the 
lattice depth changes, because for = the effect of 
the resonant and off-resonant contributions to the lower 
dressed potential are summed, whereas for — 7t/2, the 
lattice is formed solely by a resonant contribution. This is 
shown in Fig. [2^,. In addition, for = the off-resonant 
potentials shift the coupling out of resonance, changing 
the adiabatic dressed states. As a result, the admixture 
of the excited internal state in the lower dressed level 
is relatively small. As the lattices are shifted, and the 
resonant contribution dominates, the admixture of the 
excited state increases. In Fig. [2]d, we plot the admix- 
ture of the 3 Po level, averaged over one period of the 
lower dressed potential for different values of the detun- 
ing S. We note that for the detuning values we choose 
here, this value is always small. This will lead to a sig- 
nificant reduction in collisional loss rates due to 3 Po- 3 Po 
collisions when two atoms are on the same lattice site. 



D. Loading a quantum register 



In order to produce a quantum register with one atom 
in every lattice site, we begin from a spin-polarised gas of 
fermionic alkaline earth-like atoms, produced by optical 
pumping. This should be a degenerate Fermi gas so that 
the densities are sufficiently high to load a single atom 
per lattice site. Note that we choose Fermions here be- 
cause for Yb and Sr, it is the fermionic isotopes that have 
non-zero nuclear spin, and thus allow us to encode qubits 
using this degree of freedom. In the case that we have 
sufficient intensity to produce a large AC Stark shift at 
the same frequency as the final lattice, we can first load 
the gas into an off-resonance optical lattice in the 1 So 
state, and then adiabatically tune the coupling closer to 
resonance with the 3 Po state in order to load the gas 
carefully into the lower dressed potential. A high-fidelity 
quantum register can then be formed by creating a band- 
insulator state [3D], and we gain substantially over the 
case where bosons would be used for a quantum register, 
as the temperature need only be substantially smaller 
than the bandgap, and not an interaction energy for the 
band insulator to form. In addition, if a harmonic trap- 
ping potential is added to the system, most defects in the 
state will be localised near the edges of the trap [3T] . The 
resulting state can be further improved upon by apply- 
ing additional techniques, such as filtering of the state 
to improve the fidelity 32 or fault-tolerant loading of 
atoms by transfer of atoms between two internal states, 
one trapped by the lattice and the other not [33 . 
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FIG. 3: Schematic diagram of qubit readout. Qubits are 
stored in dressed states |0) and |1), which are dressed su- 
perpositions of states in the 1 So and 3 Po manifolds with a 
definite values of mi. These can be coupled via off-resonant 
Raman processes to long-lived auxilliary states |0x) and \ lx) 
in the 3 P2 manifold for the purpose of readout. In order to 
read out a particular qubit state, this state should be coupled 
to the 3 P2 manifold. It can then be detected by fluorescence 
on the cycling transition 3 P2- 3 D3. 



III. SINGLE QUBIT ADDRESSING VIA THE 
3 P 2 LEVEL 

We would like to be able to read out the state of a 
single qubit, or alternatively perform gate operations on 
a single qubit. The has been enormous recent process in 
individual addressing of sites in an optical lattice via opti- 
cal means [3"4"rl4T] . However, it would also be useful to be 
able to address individual qubits without the use of these 
techniques and the corresponding overheads in experi- 
ments. Such addressing can be achieved by coupling our 
dressed state qubit-selectively to states in the metastable 
3 P2 level, and then detecting whether the atom is indeed 
present in the 3 P2 manifold. For the purpose of readout it 
is only necessary to be able to couple one of our two qubit 
states, e.g., the |0) state (which could be represented, 
e.g., by mi = —9/2 in 87 Sr) to an auxilliary level |0x) 
in the 3 P 2 level (e.g., the | 3 P 2 , = 13/2, m F = -13/2) 
state, where F is the total angular momentum quantum 
number F and is the magnetic quantum number). 
The readout process is depicted schematically in Fig. [3] 

Because of their non-zero electron spin, states in the 
long-lived 3 P2 manifold are much more sensitive to mag- 
netic fields than the 3 Po and 1 So level, and we can use 
these shifts to make possible a spatially-dependent read- 
out of spin states by applying a magnetic gradient field, 
in a manner first mentioned in Ref. |22j . In applying such 
a field, the 3 P2 level can be significantly shifted, whilst 
the 1 So and 3 Po states are not substantially shifted, and 
thus the form of the dressed lattice potential is not sub- 
stantially changed. In particular, a gradient field of 1 
G/cm will provide an energy gradient of 4.1 MHz/cm for 



the \ 3 P2,F = 13/2, m/ = —13/2) state, or an energy 
difference of about 15 kHz between atoms in neighbour- 
ing sites for a field gradient of 100 Gauss/cm. Atoms in 
the dressed lattice can then be selectively transferred via 
a Raman process connecting off-resonantly via the 3 Si 
manifold to the 3 P2 manifold, on a timescale limited by 
the frequency shift between neighbouring sites. 

This assumes, of course, that the state in the 3 P 2 man- 
ifold to which we couple, |0x) is trapped in a lattice, 
preferably in a lattice at the same position as our qubit 
states |0) and |1). Thus, the most favourable states are 
those with a significant negative ylC-polarisability a at 
the wavelength of the clock transition, as the potential 
they experience due to the AC-Stark shift will have min- 
ima in the same places as the lower dressed state gener- 
ated by the same lattice laser. We have computed the 
polarisability from known data of the states in the 3 P2 
manifold of 87 Sr, and have found that they vary substan- 
tially due to a large tensor shift. We write the shift Ap 
from linearly polarised light as 



hA, 



-aE 2 , 



, scalar 



, tensor °" l F 



F(F + 1) 



F(2F- 1) 



E 2 



(5) 
,(6) 



where we have separated the coefficients of the scalar and 
tensor shifts 0135], and we obtain total polarisabilities at 
the clock transition frequency as shown in Fig.|4j Here we 
note that light polarised along the quantisation axis will 
give rise to a negative polarisability for the F = 13/2, 
nip = —13/2 state. This state is thus trapped by the 
same field creating the dressed lattice. We can couple 
from the mj = —9/2 states in the dressed lattice via a 
Raman process directly into the F = 13/2, mp = —13/2 
state of the 3 P2 manifold, making this state ideal for use 
as the | Ox) state in readout operations. A qubit could be 
read out by choosing the detuning of a Raman coupling 
between the | 3 P 2 , F = 13/2, m F = -13/2) state (auxil- 
iary state | Ox)) and the |— ) dressed state with mj = —9/2 
(qubit state |0)) so that it is in resonance at only one site 
as a result of a gradient field shifting the energy of the 
| 3 P 2 , F = 13/2, rap = -13/2) state. Coupling of the 
3 P2 level to a second qubit state with mj = —7/2 would 
not occur as the mp = —11/2 state is not trapped (if 
for a different species the equivalent state was trapped, 
then the large tensor shift would probably result in this 
transition being anyway out of resonance). The occu- 
pation of the 3 P2 level can then be determined by fluo- 
rescence measurements, e.g., using the cycling transition 
3 P2- 3 D3, independent of the atoms remaining in the 
and 3 Po levels. Note that the timescale for this readout 
process r rea dout is limited by the trapping frequency in 
the dressed lattice potential, T rea dout 3> 27r/w. This re- 
quirement must be fulfilled so that the atom is not cou- 
pled to excited Bloch bands of the lattice. It is also desir- 
able for this coupling to have similar trapping frequencies 
for the lattices trapping |0) and |0x), so maximising spa- 
tial overlap of the wavef unctions. Again, the F = 13/2, 
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mp = —13/2 state of the 3 P2 manifold is favourable for 
this, as the polarisability indicates that the lattice depth 
will be around 150 kHz for / ~ 3kW/cm 3 , which is a 
similar depth to that of the lattice for the dressed levels 
at the same lattice intensity (assuming that the detuning 
of the resonant coupling lasers, 5i is small). 

Note that one could equally use states with |m/| < 
— 13/2 in this process if one stores the qubit states in 
the upper dressed potential. This is disadvantageous, 
because a large detuning 6 must be chosen for the lattice 
lasers to prevent non-adiabatic loss of atoms from the 
potential [2U] . Alternatively, an additional standing wave 
at a different frequency could be added to trap states in 
3 P2 manifold via an additional AC-Stark shift. 

An alternative to using magnetic field gradients for 
addressing would be to apply a laser with spatially vary- 
ing intensity at the magic wavelength (for equal shifts of 
the 3 Po and 1 So levels. This would provide a position- 
dependent differential AC-Stark shift between the qubit 
states and the 3 P2 level, without affecting the relative 
energy of the 3 Po and 1 So levels, and thus the dressed 
lattice. 
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in the 3 P2 manifold that can be coupled to both qubit 
states. For 87 Sr, such addressing for single-qubit opera- 
tions would thus mean either using an auxiliary lattice 
to trap states from the 3 P2 manifold, or using the upper 
dressed states for qubit storage. 

Two-qubit gates can, in principle, be performed sim- 
ilarly to exisiting schemes for alkali atoms, making use 
of the spin-dependent potentials. In particular, exisit- 
ing schemes for controlled collisions can be used to pro- 
duce controlled-phase gates for atoms in neighbouring 
sites [23] . This has been implemented experimentally in a 
proof-of-principle experiment with alkali atoms [24] . but 
here we could take advantage of the 2D spin-dependent 
lattices without having to tune trapping lasers between 
fine- structure states. 

These schemes can be seen to implement controlled- 
phase gates in three steps: 

1. The spin-dependent lattices for each state are 
shifted relative to each other so that atoms at a 
chosen distance, e.g., in neighbouring lattice sites, 
will come together on the same site if and only if 
they were originally in a specific combination of 
qubit states. For example, if we write the state of 
a pair of neighbouring qubits as \q1q2), where q\ is 
the state of the first qubit and qi is the state of the 
second qubit in the pair, then atoms in the state 
1 01) are brought together, whilst |00), 1 11) and 1 10) 
remain separated (see Fig. [5| . 

2. A phase shift is generated conditioned on whether 
two atoms are on the same site or not. 

3. The atoms are returned to their initial positions. 

The phase in step two can be generated in a number of 
different methods, including via direct collisional phase 
shifts, or the use of blockade mechanisms. These differ- 
ent mechanisms are discussed in the following two sub- 
sections. 



FIG. 4: AC-Polarisabilities for 87 Sr in the 3 P 2 manifold with 
F — 13/2 at the frequency of the clock transition. 



IV. QUANTUM GATES IN SPIN-DEPENDENT 
POTENTIALS 

Single-qubit gates can be performed in one of two ways 
in this scheme. The simplest means to obtain a global 
rotation of many qubits is to directly couple the dressed 
states for two nuclear spins via a Raman process. Al- 
ternatively, different nuclear spin (qubit) states can be 
alternately coupled to auxiliary states in the 3 P2 level 
in order to provide individual addressing for single-qubit 
rotations using the techniques described in the previous 
section. Such coupling requires the use of a trapped state 



A. Phase for two-qubit gates: controlled collisions 

For alkali atoms, the phase in step 2 is generated by col- 
lisional interactions between atoms. This could be per- 
formed directly if the atoms used have a relatively large 
scattering length in the 1 So manifold (e.g., 87 Sr). For 
other species and isotopes such as 171 Yb, this could also 
be achieved using optical Feshbach resonances [J3J H3] to 
enhance the otherwise very weak collisional interaction. 
The speed of such gates is limited by the strength of 
the on-site interaction between atoms, which for a single- 
band model is limited by the trap frequency in each lat- 
tice site, lo. 

However, the existence of weak collisional interactions 
for certain isotopes also motivates us to look at other 
gate schemes, particularly using excitations to states in 
the 3 P2 manifold. 
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FIG. 5: If the qubits are trapped in a spin-dependent lattice, 
it is possible to shift the lattice for one qubit state by one 
site, so that neighbouring atoms are brought together only 
if the qubit to the left was in state |0) and the qubit to the 
right was in state |1). This can be used to aid in producing 
two-qubit quantum gates (see text for details). 



Initial State 


After Step 1 


After Step 2 


After Step 3 


10,0) 


-|0z,0a;) 


-|0x,0a;> 


|0,0) 


|o,i) 


-i\0x, 1} 


-i|0sr, 1} 


—10, 1) 


|i,o> 


-i\l,0x) 


i|l,0z> 


|1, o> 


IM) 


IM) 


— i\lx, lx) 


IM) 



TABLE I: The state of a two-qubit system after each step of 
the protocol for a blockade gate. 



the states of the two-qubit system after each step of this 
protocol are given in table [I] 

In practice, the state |0, 1) will collect a small addi- 
tional phase (f> ~ fi/A, where f2 is the coupling Rabi 
frequency and A the detuning from the excited state, 
generated by the difference between 3 Pq. 







ii) 



iii) 



B. Use of S P 2 levels 

We would again like to make use of states in the 3 P2 
manifold to which our dressed qubit states (for a fixed 
nuclear spin) can be coupled, and which are trapped in 
the same locations as our qubits. This time we will as- 
sume that we have two such auxiliary levels, \0x) and 
| la;), as depicted in Fig. [3] 



Phase for two-qubit gates: dipole blockade mechanism 



For sufficiently large onsite dipole-dipole interactions, 
which provide a energy shift between 3 P2- 3 P2 collisional 
interactions and 3 Po/ 1 So- 3 P2 corresponding to a large 
frequency shift A, we can use a dipole blockade mech- 
anism to produce a 7r phase shift, as proposed, e.g., for 
Rydberg atoms [IS]. This is illustrated in Fig. [6j and 
consists of 3 steps: 

1. Excite all |0) qubit states to an auxiliary level |0x) 
with a 7r-pulsc 

2. Couple all |1) qubit states to an auxiliary level | la;) 
with a 27r-pulse at Rabi frequency f2, assuming that 
there is no collisional interaction between the |0x) 
state and either |1) or | la;) (i.e., the pulse duration 
T is given by f2T = 2-7T. In the ideal case, if the two 
atoms are on the same site (as will happen for an 
initial state |0, 1), this step should be blocked by 
collisional interactions, which detune the coupling 
by a frequency A. 

3. Return the |0a;) state to the |0) state with a ir pulse. 

Assuming there is no coupling of the qubit state |1) to 
the auxiliary | la;) when the atom is on the same site as an 
already excited \0x) state (i.e., the blocking is perfect), 



a) 



1 So — 3 Pq 



10) 



10) 



II) 



|0) 



b) 



I to, la ) 

A—, r 
ii)/ - 



I Ok, 1)- 



OAx) 



\\x,0)- 



10,1} 



|lx,0a:) 

ii) 



|i,o) 



|l,0a:) 



FIG. 6: Schematic diagram of a blockade gate including loss, 
(a) Operations performed on the individual qubit states |0) 
and (see text for details) (b) Comparison of the opera- 
tions for initial two-qubit states |0, 1} and |1, 0} in neighbour- 
ing qubits, showing the two-qubit levels. 



2. Lossy blockade mechanism 

It was shown by C. Greene and his collaborators 
[1H1I3Z] that, in fact, two-body collisions of atoms in the 
3 P2 level lead to large inelastic loss. However, this loss 
can actually help us in producing the blockade effect, as 
large losses involving coupling to the continuum at a rate 
r from a given level can also dynamically suppress occu- 
pation of that level, as is well known from the physics of 
a two-level system. In the limit where A <C T, this would 
even produce a blockade gate based entirely on a lossy 
blockade mechanism. In this way we can turn an ap- 
parent problem into a feature of the system. Such ideas 
have also been proposed in the context of quantum sim- 
ulation with cold atoms in optical lattices, where three- 
body losses can be used to prepare interesting many-body 



states via a similar mechanism [351 |3S] 

The key characteristics of the inelastic loss processes 
that make this possible are: 

• The energy change in the inelastic collision is larger 
than the lattice depth, so that the energy carried 
away as kinetic energy is sufficient to couple the 
atoms into the continuum of motional states. 

• The length scale on which the physics of the inelas- 
tic collision takes place is smaller than the confine- 
ment length in a lattice site, so we do not expect 
the loss process to be substantially modified by the 
presence of the lattice 

• The rates for loss are large, and could reach of the 
order of T = 2ir x 20kHz for lattice densities up to 
10 16 cm- 3 . 

In the presence of loss, the basic physics of the second 
step of the protocol, as illustrated in Fig. [6] then reduces 
to a two level system, where the state with one atom in 
1 So and one in 3 P 2 playing the role of a lossless "ground" 
state and that with two atoms in 3 P2 the role of the 
lossy excited state. If we write these states as a spin-1/2 
system, the Hamiltonian reduces to 



H 



(7) 



where ct+ = \e){g\, a = \g){e\ and a z = |e)(e| - \g)(g\ 
are the usual spin operators for our two-level system with 
lossy excited state |e) and lossless "ground" state \g), ft 
is the Rabi frequency for the coupling laser, and A is the 
effective detuning from the excited state, which can be 
induced by interaction between two atoms when they are 
both in the 3 P2 manifold. Including the loss, this system 
is described by the master equation 

p = -i[H, p] - - [<T + cr~p + pa + <j- - 2a~pa + ] . (8) 

In the limit A, V 3> £1 we can describe the time evolu- 
tion of a system initially prepared in the ground state in 
perturbation theory, giving the probability that no decay 
has occurred at short times t as 



P : 



with 



n 2 



off 



4(A 2 +T 2 /4) 



r 



(9) 



(10) 



in the limit that T <C A. For our lossy blockade gate this 
is the worst-case scenario for loss events. We immedi- 
ately see that the ratio of the loss time to the gate time 
(determined by fl) is given by Ti oss /T gatc = f2/I\ This 
will give the fidelity of the lossy blockade gate. 

The blockade mechanism is illustrated in Fig. [7J where 
we plot the decay probability as a function of time i, and 
then at fixed time fit = 2n for varying T/fl. 
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FIG. 7: Loss from a two-level system prepared in the stable 
state and coupled to the lossy state computed via integration 
of eq. pj. (a) The probability that the system has under- 
gone a loss event as a function of time when prepared in the 
ground state, with A = 0, for varying values of T/Q. (b) 
The probability that a system has undergone a loss event by 
time t£l — 2ir. These values represent the gate fidelity of a 
lossy blockade gate with A = 0, or with a combined blockade 
generated by interactions and loss with A^0. 



V. OTHER TWO-QUBIT GATES 

It would also be possible to make use of exchange in- 
teractions for fermions |50) . but we will not discuss this 
in detail because it does not make specific use of the spin- 
dependent potentials, and does not, in its original form, 
take specific advantage of the properties of alkaline earth 
atoms. 

Another possibility is the direct use of Rydberg gates 
[45] . which have been recently demonstrated for trapped 
alkali atoms [5"T1 152] . The separate hierarchy of Ryd- 
berg states for the singlet and triplet manifolds could 
give advantages for Rydberg excitations in alkaline earth 
atoms, especially facilitating easier state-dependent ex- 
citation. These could also be performed together with 
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gradient addressing, exciting the Rydberg state from the 
3 P2 manifold. 



VI. DECOHERENCE/LOSS MECHANISMS 

There will be a number of possible sources of decohcr- 
ence within this setup, all of which should be controllable 
in the experiment. These include magnetic field fluctu- 
ations, decoherence due to frequency noise of the lasers, 
spontaneous emissions, and collisional losses. We briefly 
summarise the role of these key sources of decoherence 
below. 



A. Magnetic field fluctuations 

Magnetic field fluctuations will probably still consti- 
tute the largest source of decoherence, however, this is 
reduced by almost 3 orders of magnitude compared with 
qubits encoded on an electron spin. Decoherence from lo- 
cal fluctuations in the magnetic field will contribute both 
due to direct shifts of the energy of the qubit states, and 
from the modification to the dressed potentials due to 
the differential shift between the 1 Sq and 3 P levels. 



as for fli ^> Si AV « ASf /CI 2 , this is also suppressed by 
an extra factor of AS /CI. In addition, other sources of 
heating, such as intensity noise on the lasers creating the 
lattice, or shaking of the lattice potential (due, e.g., to 
vibrating optical components) will have a similar effect 
[531 [54]. 



C. Spontaneous emissions 

Qubits can decohere or be destroyed (the atoms lost 
from the lattice) by spontaneous emission events. These 
can come from two sources: the finite lifetime of the 3 Po 
state, and off-resonant coupling to states with a short life- 
time induced by the lattice lasers (storage) or coupling 
lasers (during gate and readout operations). However, 
the lifetime of 3 Po is many seconds, and this source of 
atom loss can be further suppressed by using the resonant 
lattices only for spin-dependent transfer, and storing the 
atoms at other times in the state (see previous sub- 
section) . 



D. Collisional losses from 3 Pq 



B. Stability of the trapping lasers 

A finite laser linewidth for the dressing laser creating 
the potentials will give rise to fluctuations ASi in the 
detuning Si, and therefore the energy of atoms trapped 
in the dressed potential. However, these fluctuations will 
lead to the same fluctuation in lattice depth for the two 
qubits, ASq — AS i . The resulting ground state energy 
will shift by different amounts, as the lattice periods are 
different. However, if the corresponding wavelengths are 
A and (1 + e)X, then the difference in trap frequencies for 
the two qubit states, A u is given in terms of the depth 
fluctuations AV^ by 



^ = Jav*^ 

2 V 2mA 2 



' AV4tt 2 H 2 
2mX 2 (l + e) 2 



-e\ SV 



4TT 2 h 2 



2mA 2 
(11) 

Thus, as e ~ 10 8 , this decoherence mechanism will be 
strongly suppressed, and for laser linewidths of the or- 
der of tens of Hz, dephasing times can be many min- 
utes. On the other hand, the resulting noise AV on the 
depth of the lattice could give rise to heating of the par- 
ticles to higher oscillator levels, if appropriate frequency 
components are present in the noise in order to drive 
these coupling. s This would lead to imperfect couplings 
for gate and readout operations. Such heating rates 
can be estimated [531 154] as giving an energy increase 
(E) = r heat {E) with rate r heat = it 2 lo 2 S e {2uj) /2, where 
S e (2u) is the one-sided power spectrum of the trap am- 
plitude noise at twice the trap frequency ui. In our case, 



Measurements of collisional losses between atoms in 
the 3 Po manifold are currently underway in several 
groups, in order to determine what the collisional life- 
time is when two atoms are present in these states at 
the typical lattice densities that will be encountered here 
(ca. 10 14 cm~~ 3 -10 15 cm~ 3 onsite). Effects of these losses 
have been observed recently, e.g., in samples of Stron- 
tium atoms confined in ID tubes [55] - However in our 
case, during storage, readout, and single-qubit opera- 
tions, the atoms are anyway isolated by the lattice, and 
two atoms will not collide. Thus, the only time that two 
atoms with components of states from the 3 Po manifold 
are present on the same site is during two-qubit gate op- 
erations. If these take place on a timescale ca. 1ms, 
then we would require collisional stability of our atoms 
for timescales longer than 100ms in order to achieve gate 
fidelities larger than 99% if both atoms were in the 3 Po 
manifold. However, as shown in Sec. |II C| the combina- 
tion of resonant and off-resonant lattices mean that the 
amplitude for atoms to be in the 3 Po manifold is small for 
all stages of operation except during transport of atoms. 
Gate schemes can be made more immune to these losses 
by using larger intensity trapping lasers, and thus intro- 
ducing a larger component from the off-resonant lattice 
(see Sec. II C ). This will ensure that when the lattices for 



the two qubit states overlap that the dressed states are 
dominated by off-resonant lattices for 1 So, and that the 
admixture of the 3 Po state is small. If the probability to 
find a single atom on a given site in the 3 Po manifold is 
e 3 for each of the qubit states, then the onsite loss rate 
will be suppressed by a factor ~ e 2 . 
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VII. SUMMARY AND OUTLOOK 

In summary, the quantum computing scheme we pre- 
sented based on nuclear-spin-depenent lattices with near- 
resonant coupling on the clock transition for alkaline- 
earth(-like) atoms has several advantages over schemes 
with alkali atoms. The use of nuclear spins for qubit 
storage makes this scheme relatively robust against de- 
coherence due to magnetic field fluctuations, and cou- 
pling to the 3 P2 manifold provides high-resolution indi- 
vidual qubit addressing with a magnetic gradient field. 
There are also possibilities here to perform gates based 
on transfer of states to long-lived metastable excited lev- 
els (e.g., 3 P2), including the new mechanism of lossy 
blockade gates. In comparison with a scheme presented 
previously using electronic-state-dependent lattices, this 
scheme does not require lasers that independently trap 
the 1 So and 3 Po manifolds. This method is more sensitive 
to collisional losses between two atoms in the 3 Po mani- 
fold, although this only affects short periods of time dur- 
ing the gate operations. While we have focused here on 
gates based on state-dependent lattices, other schemes, 
including Rydberg gates will benefit from the unique 
properties of alkaline-earth-like atoms. In particular, 
state-selective excitation to a Rydberg state would be 
simplified, e.g., by exciting one nuclear spin state to the 
3 P 2 manifold first. 

The key experimental requirements for implemention 
of these methods are: (i) Large, stable magnetic fields 
(to provide the differential Zeeman shifts allowing spin- 
dependent lattices for different nuclear spin states; (ii) 
High-intensity stable laser on the clock transition (to pro- 



vide a deep optical lattice whilst avoiding decoherence 
due to noise on the detuning 8q.); and (iii) Control over 
magnetic field gradients (to allow for either large paral- 
lel operations or individual addressing with qubits oper- 
ations involving coupling to 3 P 2 (although single-qubit 
gates could also be done directly in parallel, and the use 
of 3 P 2 is only necessary in two-qubit gates in the case 
that the scattering lengths for the clock states are not 
sufficiently large). 

Nuclear-spin-dependent lattices also have immedi- 
ate possible application for quantum simulation with 
alkaline-earth-metal atoms. In particular, the depen- 
dence on the nuclear spin state could be used to break 
the degeneracy in models with SU(N) symmetry (251 RZo] , 
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